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'  ABSTRACT 

tf 

There  are  many  physical  sources  that 
generate  recurring  pulses.  Rotating  fans  and 
propellers  are  two  examples.  The  objective  of 
this  paper  is  to  present  a  signal  processing 
methodology  to  extract  recurring  pulses  from 
data.  The  method  relies  on  obtaining  discrete 
frequency  components  over  a  band  or  discrete 
temporal  components  over  a  time  interval.  It  is 
shown  that  the  results  for  the  time  and  frequency 
domains  are  duals  of  each  other.  These  ideas  are 
then  generalized  to  the  spatial  domain  by 
including  an  array  of  sensors.  It  is  shown  that 
the  time-spatial  domain  and  the 
frequency-wavenumber  domain  are  also  duals  of 
each  other. 

A. 

INTRODUCTION 

The  physical  principle  of  extracting 
recurring  pulses  from  data  relies  on  obtaining 
discrete  discernable  components  either  in  the 
frequency  or  time  domain.  Data  from  rotating 
machinery  usually  produce  discrete  harmonically 
related  frequency  components.  On  the  other  hand, 
discrete  recurring  temporal  pulses  are  produced 
by  other  sources.  For  example,  helicopter- 
radiated  noise  [1]. 

We  shall  first  consider  a  signal  processing 
methodology  for  a  single  sensor,  to  extract 
recurring  pulses  or  harmonically  related 
frequency  components.  The  results  for  these  two 
cases  are  related  through  a  time-frequency 
duality  principle.  In  addition,  the  extracted 
components  are  equivalent  to  the  maximum 
likelihood  estimate.  On  the  other  hand,  for  an 
array  of  sensors  the  signal  processing  meti.od  for 
these  two  cases  are  related  through  a 
t ime-spati al -frequency-wavenumber  dual  1 ty 
principle. 

In  the  paper  the  data  are  assumed  to  be 
filtered  to  the  desired  passband  and 
appropriately  A/D  converted. 


The  discrete  Fourier  transform  is  defined  as 
follows. 


X(K)  .  (1/Vn)  jp  x ( i )  exp(-j2»Ki/N) 


and  its  inverse 

N-l 

x(i )  -  (llyfn)  £  X(K)  exp(j2»Ki/N) 

K-o 

For  an  interval  of  time  consisting  of  N 
samples  let  the  data  be  represented  by  a 
succession  of  pulses  recurring  every  r  seconds. 
This  can  be  written  as  follows. 


x ( i )  -  ^  «[i  -  (ig  ♦  kr)] 


Where  «(  )  is  the  Kronecker  delta  function. 
Here,  ig  represents  the  sample  on  which  the 
first  of  w  pulses  starts. 

Notice  that  the  spectrum  of  the  above 
process  is  given  by 


X(K)  -  {w lyfS)  exp[-J2.(i  ♦  )  K/N] 


I  sin(«wlWN)  | 

[w  s i n ( wK T /N ) J  * 

The  location  of  the  peaks  of  this  spectrum 
are  found  from  the  solutions  to  the  equation 


sin(wlWN)  -  0 


which  are  given  by,  K  •  (N/r)l0,  for 
o  <_  K  ^  N-l,  and  Ig  «  0,1,2,  ... 


1 


w: 


2 


Therefore ,  r  controls  both  the  number  and 
locations  of  the  discrete  components.  As  w 
the  spectrum  of  the  process  will  approach  a 
purely  line  spectrum. 

This  example  reveals  the  principle  on  which 
the  signal  processing  method  is  based.  Namely, 
all  the  information  associated  with  recurring 
pulses  is  contained  at  discrete  frequencies.  The 
dual  of  this  result  is  that  all  the  information 
associated  with  the  recurring  pulses  is  contained 
at  discrete  time  locations. 


MATRIX  METHOD 

The  discrete  Fourier  transform  and  its 
inverse  can  be  conveniently  discussed  in  terms  of 
matrices  and  vectors. 


x.  -  (*(0),  x(l),  ....  x(N-l) )T 


be  a  vector  of  time  domain  data  samples.  Here  T 
represents  transpose. 

Similarly,  the  frequency  domain  data  is 
given  by  the  vector 


X  -  (X(0),  X ( 1 ) . X(N-1))t 


where  the  two  vectors  are  related  through  the 
transformation 


X  -  F  x 


and  F  is  a  (NxN)  transform  matrix  representing 
the  OFT  operation  as  follows: 


We  will  now  define  a  diagonal  (NxN)  operator 
matrix  Q  which  is  applied  in  the  frequency  domain. 

The  extracted  time  domain  data  vector  x£, 
after  the  matrix  Q  is  applied  in  the  frequency 
domain  is  given  by  the  expression 


X£  -  F-lQF  x 


Notice  that  the  time  domain  transformation 
F~1qF  is  real.  If  we  subtract  X£  from  x  we 
obtain  the  resultant  data  vector,  xg. 


xg  »  x  -  F-lQF  x. 


If  Q«I  (identity  matrix),  then,  xg  *  0  (zero 
vector). 

Whereas,  if  Q  »  0  (zero  matrix),  then  xg  «  x. 

Therefore  we  can  choose  Q  to  obtain  the  desired 
result. 

In  applications  the  elements  of  the  diagonal 
matrix  Q  are  chosen  to  represent  the  locations 
where  the  discrete  components  occur  in  the 
frequency  domain. 

For  example,  let  N-1024,  and  let  there  be 
four  exact  solutions  located  at  K-0,  256,  512, 
and  768. 


Then  Q,  -  diag  [q^-;  i,j«0,l .  1023] 


where. 


1,  at  i-j  -  0,  256,  512,  and  768 


0,  otherwise 


(1/N)1'2 


1  1 
1  exp(-j2*/N) 

1  exp(-j2*2/N) 


exp(-j2*2/N) 

exp(-j2»4/N) 


exp(-j2*(N-l)/N) 

exp(-J2*2(N-l)/N) 


1  exp(-j2*(N-2)/N)  exp(-j2»2(N-2)/N  . 
j_l  exp(-j2*(N-l)/N)  exp(-j2»2(N-l)/N)  . 

and  its  inverse  exists  and  is  denoted  by 


F'1  »  (l/N)l/2 


1  1  1 
1  exp( j2»/N)  exp( J2*2/N) 
1  exp( j2»2/N)  exp{ J2*4/N) 


exp(-j2t(N-2)(N-l)/N) 
exp(-j2«(N~l ){H-1 )/N) 


exp( j2#(N-l )/N) 
exp(j2*2(N-l)/N) 


l  exp(j2*(N-2)/N)  exp( j2*2(N-2) 7N) 
1  exp( J2*(N-1)/N)  exo( j2»2(N-l )/N) 


exp{ j2*(N-2)(N-l )/N) 
exp( J2»(N-1)(N-1)/N) 


■i  i it-; * ;  ■  :• 


Therefore,  for  this  example,  the  extraction 
method  is  nearly  perfect  even  with  additive 
noise.  In  practice  the  matrix  method  will  be 
limited  by  the  ability  to  choose  the  appropriate 
operator  matrix  Q. 

We  shall  now  compare  this  method  witn  the 
maximum  likelihood  (M.)  estimate  [2].  Let, 


2  m  x  ♦  n 


be  the  data.  Here  rt  is  i.i.d.  Gaussian  noise. 

We  want  to  est1mate~the  deterministic  waveform  x 
at  each  sample  point.  These  sample  points  will 
be  denoted  by  a.  The.  ML  estimate  of  a  is  given 
by  the  procedure.  Choose  the  value  oT  a,  say  $, 
that  maximizes  f(z\i),  where,  here,  f(  T  is  a 
multidimensional  Eaussian  probability  density 
function.  Since,  the  waveform  is  a  recurring 
pulse  of  duration  h,  the  estimate  reduces  to  the 
expression 

w-1 

*1  ■  U/")  £  x(ig+i+hk) 
k-0 


Xg  -  (4/N)  £  X  a  X 


So  that,  tne  extracted  data  vector  is  equal 
to  the  original  data  vector.  The  reason  for  this 
is  that,  in  this  example,  all  the  information  was 
located  at  discrete  frequencies.  Now,  if  we  add 
noise  of  the  form  n  then 


Z  m  X  *  £ 

and 

2£  -  (4/N)  Ex*  (4/N)  E  n 
-  x  ♦  (4/N)  E  n 


which  reduces  to,  Z£  -  x.  Since  the  noise 
is  random,  the  noise  vector  will  not  match  the 
transformation  E  and  consequently  the  noise  will 
be  effectively  averaged  out. 

The  resultant  data  vector  is  obtained  from 
the  expression. 


XR  «  x_  -  (4/N)  E  x  -  0. 


Adding  noise  again  we  obtain, 


Zfl  «  n  -  (4/N)  £  n  -  n  . 


for  i-0,l,2...,h-l.  So,  each  sample  of  tne  pulse 
is  averaged  over  all  respective  samples  of  the 
successive  pulses.  This  is  exactly  what  tne 
transformation  F-lQF  performs. 

For  a  sinusoid,  which  is  periodic,  tne  same 
type  of  transformation,  but  dependent  upon 
frequency,  will  evolve  from  tne  extraction 
procedure.  Here,  however,  two  components  are 
associated  with  each  sinusoid  [3]. 


For  example,  let  the  frequency  of  a  sinusoid 
be  k>256.  Then  the  operator  matrix  Q  will  have 
two  nonzero  components  located  at  q,j  (i-j-256 
and  768).  The  transformation  will  reduce  to  the 
following. 


It  is  clear  from  this  result  that  the  exact 
sinusoid,  including  phase,  will  be  extracted. 

Now  consider  a  frequency  domain 
transformation.  This  procedure  is  a  dual  of  the 
time  domain  transformation. 


The  extracted  frequency  domain  data  vector, 
Xg  is  given  by  the  expression 

X£  .  FQtF_1X 


SPATIAL  DOMAIN 


These  results  can  be  generalized  for  an 
array  of  sensors.  The  two-dimensional  discrete 
Fourier  transform  is  given  by  the  expression 


Where  Qt  is  a  diagonal  operator  matrix  applied 
in  the  time  domain.  Now,  FQtF"*  will  be 
complex  in  general.  But  the  transformation 
performs  a  similar  operation  as  its  time  domain 
dual.  Namely,  it  averages  the  appropriate 
frequencies  to  produce  the  extracted  data 
vector.  This  is  also  equivalent,  under  the 
stated  conditions,  to  the  ML  estimate. 

The  resultant  frequency  domain  data  vector, 
Xg,  is  given  by  the  expression. 


Xg  -  X  -  FQtF-lx 


Transforming  this  result  into  the  time 
domain  we  obtain  the  result. 


*R  -  [I-Qtl  *. 


The  transformation  [I-Qtl  is  equivalent  to 
the  ideal-nonlinearity  discussed  in  [3,1].  This 
means  that  if  we  can  identify  the  recurring 
pulses,  or  if  we  have  a  priori  knowledge  of  them, 
in  the  time  domain  and  apply  the  transformation 
[I-Qtl  we  will  have  subtracted  the  frequency 
domain  ML  estimate  from  X. 

Suppose  the  nonzero  components  of  Qt  are 
located  at  the  sample  points,  0,  256,  512  and 
768.  Then  the  transformation  will  also  reduce  to 
FQtF-1  -  (4/N)E. 

An  important  application  of  this  method  is 
obtaining  the  power  spectrum  of  a  signal  when 
recurring  impulsive  interference  is  present. 

The  resultant  data  vector  will  demonstrate  the 
ability  to  remove  this  Interference. 

From  the  above  result  we  obtain  the  power 
spectrum  from  the  the  expression 


M-l  N-l 

X(L,K)  .  (lAjM) (1/v'N)  £  £ 

J«0  i»0 

exp[-j2«(JL/M  -  iK/N) ] 

where,  J-0,1,...,M-1  represent  the  spatial 
sensors,  i.e.,  x(0,i)  is  the  first  sensor  at  the 
i-th  time  sample.  If,  N-l,  x(J,0)  and  X(L)  will 
have  the  same  relationship  as  x(i)  and  X(K)  did 
before.  But,  now,  J  represents  spatial  position 
and  L  spatial  frequency. 

It  is  convenient  to  write  this  equation  in 
matrix  form. 


JSFW  -  fFW  sFW  ITS 

where,  Xpy  1s  a  (MNxl )  vector  and  Fpw  is  a 
block  diagonal  (MNxMN)  matrix.  Each  block  is  the 
(NxN)  matrix  F.  The  other  matrix  Spy  is  a 
(MX  x  MN)  matrix  of  the  form. 


I  I 

I  Ie-J2»/M 


‘le-j2ir(M-l)/M 


(1/vM)  . 


|_i  Ie-j2«(M-l)/M  .  .  . Ie-j2w(M-l ) (M-l )/MJ 

where  I  is  a  (NxN)  identity  matrix,  and  xj$  is 
a  (MNxl)  vector. 

We  now  define  a  (MNxMN)  diagonal  operator 
matrix  Qpy,  so  that 


5TS  *  SFW  FFW  QFW  FFW  SFW  iTS 


*T  —1 

XR '  XR  -  X  X  -  X  FQtF  *X 


Therefore,  the  resultant  power  spectrum  is 
obtained  by  subtracting  the  power  spectrum  of 
Xp  from  the  power  spectrum  of  the  original 
<Jata.  An  error  analysis  of  this  procedure  was 
discussed  in  [4]. 


ts  the  extracted  data  vector  of  the  time-space 
domain. 

The  dual  of  this  result  is  the  following. 


^FW  “  FFW  SFW  Ots  -F«  ’ 


Where  is  a  (MNxMN)  diagonal  operator  matrix 
applied  in  the  time-space  domain. 
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SUMMARY 

Since  there  are  many  physical  sources  that 
produce  discrete  spectra  or  discrete  temporal 
components,  a  signal  processing  methodology  Mas 
developed  to  take  advantage  of  these  phenomena. 

If  a  source  produced  discrete  frequencies  then 
all  of  the  information  associated  Mith  the  source 
is  contained  at  those  discrete  frequencies.  By 
extracting  only  those  frequencies  a  time  domain 
enhancement  Mas  possible.  By  considering  the 
transformation  of  the  method  it  Mas  revealed  that 
the  extracted  data  is  equivalent  to  a  maximum 
likelihood  estimate  under  the  assumption  of 
i.i.d.  Gaussian  noise.  A  dual  relationship  Mas 
shown  to  hold  for  discrete  temporal  components. 
These  results  were  then  generalized  for  an  array 
of  sensors.  In  this  case  the  results  revealed  a 
time-space-frequency-wavenumber  duality. 
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